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Abstract 

All existing positive results on two neutrino double beta decay in different 
nuclei were analyzed. Using the procedure recommended by the Particle Data 
Group, weighted average values for half-lives of 48 Ca, ' 6 Ge, 82 Se, 96 Zr, 100 Mo, 
100 Mo - 100 Ru (Of), 116 Cd, 130 Te, 136 Xe, 150 Nd, 150 Nd - 150 Sm (Of) and 238 U 
were obtained. Existing geochemical data were analyzed and recommended 
values for half-lives of 128 Te, and 130 Ba are proposed. Given the measured 
half-life values, nuclear matrix elements were calculated using latest (more 
reliable and precise) values for phase space factor. Finally, previous results 
(PRC 81 (2010) 035501) were up-dated and results for 136 Xe were added. 

Keywords: Double beta decay, Nuclear matrix elements, 48 Ca, 76 Ge, 82 Se, 
96 Zr, 100 Mo, 100 Mo - 100 Ru (Of), 116 Cd, 128 Te, 130 Te, 136 Xe, 150 Nd, 150 Nd - 
150 Sm (Of), 238 U and 130 Ba 


1. Introduction 

At present, the two neutrino double beta (2 ufdfd) decay process has been 
detected in a total of 11 different nuclei. In 100 Mo and 150 Nd, this type 
of decay was also detected for the transition to the 0 + excited state of the 
daughter nucleus. For the case of the 130 Ba nucleus, evidence for the two 
neutrino double electron capture process was observed via a geochemical 
experiments. All of these results were obtained in a few tens of geochemical 
experiments and more than forty direct (counting) experiments as well as 
and in one radiochemical experiment. In direct experiments, for some nuclei 
there are as many as eight independent positive results (e.g., 100 Mo). In 
some experiments, the statistical error does not always play the primary 
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role in overall half-life uncertainties. For example, the NEMO-3 experiment 
with 100 Mo has currently detected more than 219,000 2 1//3(3 events |l|, which 
results in a value for the statistical error of ~ 0.2% . At the same time, 
the systematic error for many experiments on 2 v/3/3 decay remains quite 
high (~ 10 — 30%) and very often cannot be determined reliably. As a 
consequence, it is frequently difficult for the “user” to select the “best” half- 
life value among the results. Using an averaging procedure, one can produce 
the most reliable and accurate half-life values for each isotope. 

Why are accurate half-life periods necessary? The most important moti¬ 
vations are the following: 

1) Nuclear spectroscopy. Now we know that some isotopes which were earlier 
considered to be stable are not, and decay via the double beta decay pro¬ 
cesses with a half-life period of ~ 10 18 — 10 24 yr are observed. The values 
which are presented here should be introduced into the isotope tables. 

2) Nuclear matrix elements (NME). First, it gives the possibility to im¬ 
prove the quality of NME calculations for two neutrino double beta decay, so 
one can directly compare experimental and calculated values. For example, 
so-called '' g a (axial-vector coupling constant) quenching” problem could be 
solved by comparison of exact experimental values of NMEs and results of 
theoretical calculations (see discussions in Ref. 0, S 0,0. Second, it gives 
the possibility to improve the quality of NME calculations for neutrinoless 
double beta decay. The accurate half-life values for 2 v/3{3 decay are used to 
adjust the most relevant parameter of the quasiparticle random-phase ap¬ 
proximation (QRPA) model, the strength of the particle-particle interaction 


9PP . — , 

3) Research on the single state dominance (SSD) mechanism 11(1 lljand 
a check of the ’’bosonic” component of the neutrino hypothesis mm] is 
possible. 

In the present work, an analysis of all “positive” experimental results 
has been performed, and averaged or recommended values for isotopes are 
presented. 

The first time that this work was done was in 2001, and the results were 
presented at the International Workshop on the calculation of double beta 
decay nuclear matrix elements, MEDEX’01 [l^]. Then revised half-life values 
were presented at MEDEX’05 and MEDEX’09 and published in Ref. 
and 111 


17], respectively. In the present paper, new positive results obtained 


15 


since end of 2009 and to the end of 2014 have been added and analyzed. 
The main differences from the previous analysis 0 are the following: 
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1) The new experimental data obtained after the publication of Ref 


are included in the analysis: 76 Ge |l8j|, 100 Mo 
116 Cd 0, 130 Te 0, 150 Nd - 150 Sm (Of) E 



100 Mo 

and 130 Ba0 


100 


Ru(0f) 



are analyzed. This decay 


2) Positive results obtained for 136 Xe 
was detected for the fist time in 2011 0]. 

3) To calculate NMEs new phase space factor values (G 2|/ ) from Ref. 
0, 30] and 0, 32] are used. 

4) Considering possible changes of axial vector coupling constant qa (pos¬ 
sible quenching effect in nuclear medium) so-called effective NMEs are cal¬ 
culated, | M 2 I |= g\- | ( m e c 2 
matrix elements | (m e c 2 • 


M, 


2u 


in Ref. 17 the dimensionless nuclear 


were calculated for qa = 1.254). 


2. Present experimental data 

Experimental results on 2decay in different nuclei are presented in 
Table 1. For direct experiments, the number of events and the signal-to- 
background ratio are presented. 


3. Data analysis 

To obtain an average of the ensemble of available data, a standard weighted 
least squares procedure, as recommended by the Particle Data Group 0], 
was used. The weighted average and the corresponding error were calculated, 

as follows: _ 

x ± 8x = ^WiXi/^Wi ± (J^Wi)" 1/2 , (1) 

where Wi = 1 /(Ac;) 2 . Here, x t and 5xi are the value and error reported by 
the i-th experiment, and the summations run over the N experiments. 

The next step is to calculate y 2 = w i( x ~ x i ) 2 and compare it with 
N - 1, which is the expectation value of y 2 if the measurements are from a 
Gaussian distribution. If y 2 /(lV — 1) is less than or equal to 1 and there 
are no known problems with the data, then one accepts the results to be 
sound. If y 2 /(N — 1) is very large (>> 1), one chooses not to use the 
average. Alternatively, one may quote the calculated average, while making 
an educated guess of the error, using a conservative estimate designed to 
take into account known problems with the data. Finally, if y 2 /(iV — 1) is 
larger than 1, but not greatly so, it is still best to use the average data, but 
to increase the quoted error, 8x in Equation 1, by a factor of S defined by 

s = {xV(N - l)p 2 . ( 2 ) 
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Table 1: Present, positive 2vj3 /? decay results. Here N is the number of useful events, TR 2 
is a half-life and S/B is the signal-to-background ratio. 

For E 2 e >1.2 MeV; after correction (see text); c) for the SSD mechanism; ^ for 
i? 2 e > 1-5 MeV; e ) in both peaks. 


Nucleus 

N 

T 1 / 2 , yr 

S/B 

Ref., year 

48 Ca 

~ 100 

5 

116 

[4.3/ 2 ; 4 (staf) ± lA(syst)} ■ 10 19 
4.2+ 8 ; 8 • 10 19 

[4A+Q'l(stat) ± OA(syst) ■ 10 19 

Average value: 4.4l° ® • 10 19 

1/5 

5/0 

6.8 

[22], 1996 
[24], 2000 
[35], 2011 

76 Ge 

~ 4000 

758 

- 330 
132 

~ 3000 

~ 80000 

7030 

(0.9 ±0.1) • 10 21 
lAto 6 3 • 10 21 

0.921o'q4 • 10 21 

1.27+oii6 • 10 21 
(1.45 ±0.15) • 10 21 
[1.74 ± 0.01(stat))(o'}g(st/st)] • 10 21 
1.841 9 ;} 4 -'lO 21 

Average value: 1. 65^12 ' 10 21 

~ 1/8 
~ 1/6 
~ 1.2 

~ 1.4 

~ 1.5 

~ 1.5 

[36] , 1990 

[37] , 1991 
[22], 1991 
[22], 1994 

[40] , 1999 

[41] , 2003 
[12], 2013 

82 Se 

89.6 

149.1 

2750 

1-O8lo!o6 • 1Q2 ° 

[0.83 ± OAO(stat) ± 0.07(sj/s*)] • 10 2 ° 
[0.96 ± 0.03(stat) ± OA(syst)] • 10 2 ° 
(1.3 ± 0.05) • 10 2 ° (geochem.) 

Average value: (0.92 ± 0.07) • 10 2 ° 

~ 8 

2.3 

4 

[42] , 1992 

[43] , 1998 
[1], 2005 

[44] , 1986 

96 Zr 

26.7 

453 

[2.l/ 9 ; 8 (sfat) ± 0.2 (syst)} • 10 19 
[2.35 ± OAA(stat) ± 0A6(syst)] • 10 19 
(3.9 ± 0.9) • 10 19 (geochem.) 

(0.94 ± 0.32) • 10 19 (geochem.) 

Average value: (2.3 ± 0.2) • 10 19 

1.9 a) 

1 

[45] , 1999 

[46] , 2010 

[47] , 1993 

[48] , 2001 
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Table 1: continued. 


100 Mo 

~ 500 

11.5±l;g • 10 18 

1/7 

[49], 1991 


67 

11.6±g;| • 10 18 

7 

[50], 1991 


1433 

[7.3 ± 0.35 (stat) ± 0.8(syst)] ■ 10 186 ) 

3 

[51], 1995 


175 

7.6■ 10 18 

1/2 

[52], 1997 


377 

[6-75 toil(stat) ± 0.68(st/st)] • 10 18 

10 

[53], 1997 


800 

[7.2 ± 1.1 (stat) ± 1.8 (syst)] ■ 10 18 

1/9 

[54], 2001 


219000 

[7.11 ± 0.02 (stat) ± OM(syst)] • 10 18c > 

40 

yy, 2005 


~ 350 

[7.15 ± 0.37 (stat) ± 0.66 (syst)] ■ 10 18 

~ 5 d ) 

[19], 2014 



(2.1 ± 0.3) • 10 18 (geochem.) 


[55], 2004 



Average value: (7.1 ± 0.4) • 10 18 



100 Mo - 

133 e) 

6.ll]; 8 • 10 2U 

1/7 

[56], 1995 

100 Ru (0^) 

153 e ) 

[9.3± 1.4(st/s£)] • 10 20 

1/4 

[57], 1999 


19.5 

[5.9l};i(siai) ± 0.6(sj/st)] • 10 20 

~ 8 

[58], 2001 


35.5 

[5.5lJ;g(siai) ± 0.3 (syst)] ■ 10 20 

~ 8 

[59], 2009 


37.5 

[5.7jlJ;g(stat) ± 0.8(s?/st)] • 10 20 

~ 3 

[60], 2007 


597 e ) 

[6.9lo;g(stat) ± 0.7(syst)] ■ 10 20 

~ 1/10 

[20], 2010 


239 e ) 

[7.5 ± 0.6 (stat) ± 0.6 (syst)] ■ 10 20 

Average value: 6.7lo ® ■ 10 20 

2 

[21], 2014 

116 Cd 

~ 180 

2.6^| • 10 19 

~ 1/4 

[61], 1995 


174.6 

[2.9 ± 0.3 (stat) ± 0.2 (syst)] • 10 19fe ) 

3 

[62], 1996 


9850 

[2.9 ± 0.06 (st at) (syst)] ■ 10 19 

~ 3 

[63], 2003 


7000 

[2.88 ± OM(stat) ± 0.16 (syst)] • 10 19c) 

10 

[64], 2011 


34927 

[2.80 ± 0.05(stat) ± 0.4(st/st)] • 10 19 

Average value: (2.87 ± 0.13) ■ 10 19 

2 

[22], 2014 

128 Te 


~ 2.2 • 10 24 (geochem.) 


[65], 1991 



(7.7 ± 0.4) • 10 24 (geochem.) 


[66], 1993 



(2.41 ± 0.39) • 10 24 (geochem.) 


[67], 2008 



(2.3 ± 0.3) • 10 24 (geochem.) 


[68], 2008 



Recommended value: (2.0 ± 0.3) ■ 10 24 
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Table 1: continued 2. 


130 Te 

260 

236 

[6.1 ± lA(stat)+l' 9 5 (syst)] ■ 10 20 
[7.0 ± 0.9 (stat) ± 1.1 (si/si)] • 10 2 ° 

— 8 • 10 20 (geochem.) 

(27 ± 1) • 10 2 ° (geochem.) 

(9.0 ± 1.4) • 10 20 (geochem.) 

(8.0 ± 1.1) ■ 10 20 (geochem.) 

Average value: (6.9 ± 1.3) • 10 20 

1/8 

1/3 

[69], 2003 
[23], 2011 

[65] , 1991 

[66] , 1993 

[67] , 2008 

[68] , 2008 

136 Xe 

- 50000 

19042 

[2.30 ± 0.02(stat) ± 0.12(syst)} • 10 21 
[2.165 ± 0.016 (stat) ± 0.059(s?/st)] • 10 21 

Average value: (2.19 ± 0.06) • 10 21 

-10 

-10 

[26], 2012 
[22], 2014 

150 Nd 

23 

414 

2018 

[18.8l|;g(siai) ± 1.9(st/s£)] • 10 18 
[6.75lJJ;42 (stat) ± 0.68 (syst)] ■ 10 18 
[9.11 tiH(stat) ± 0.63(syst)\ ■ 10 18 

Average value: (8.2 ± 0.9) • 10 18 

1.8 

6 

2.8 

[70], 1995 
[53], 1997 
[21], 2008 

150 Nd - 

150 Sm (0/) 

177.5 d ) 

21.6 

[1-33 t° 0 f 3 (stat)toil(syst)] ■ 10 20 
[1.07 tof 5 ( s tdt) ± +0.07(si/st)] • 10 2 ° 

Average value: 1.21q 2 ' 10 20 

1/5 
- 1.2 

[72], 2009 
[24], 2014 

238-q 


(2.0 ± 0.6) • 10 21 (radiochem.) 


[ 22 ], 1991 

130 Ba 

ECEC(2z/) 


2 .1 ^o'8 ‘ 10 21 (geochem.) 

(2.2 ± 0.5) ■ 10 21 (geochem.) 

(0.60 ± 0.11) • 10 21 (geochem.) 

Recommended value: ~ 10 21 


[ 24 ] , 1996 
[75], 2001 

[25] , 2009 
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For averages, the statistical and systematic errors are treated in quadrature 
and used as a combined error Sx t . In some cases, only the results obtained 
with high enough signal-to-background ratio were used. 


3.1. 48 (7a 

There are three independent experiments in which 2 v/3/3 decay of 48 Ca 


was observed 33j, [34|, 135]. The results are in good agreement. The weighted 


average value is: 


Ti/ 2 — 4.4 


+0.6 

- 0.5 


• 10 19 yr. 


3.2. 76 Ge 

Considering the results of six experiments, a few additional comments are 
necessary, as follows: 

1) We use here final result of the Heidelberg-Moscow Collaboration, Xi/ 2 = 
[1.74 ± 0.01(stat)^o;} 8 (s2/st)] • 10 21 yr 


2) In Ref. [38j], the value T\/ 2 = 0.92/q;q4-10 21 


yr was presented. However, 
after a more careful analysis, this result has been changed ( 39 J. In Ref. 39 


a few values for half-life using different analysis methods were obtained. I 
use here the value obtained by fit the data using y 2 model, which take into 
account shape of the spectrum (see Table 1). Unfortunately systematic error 
was not discussed and taken into account in this paper. This is why during 
my analysis I added typical systematic error for such sort of experiments 
(±10%). So, finally, I use Ti/ 2 = [1.27lQ 21 6 (sfat) ± 0.13(s?/sf)] • 10 21 yr as a 
result of Ref. Q. 

3) The results presented in Ref. |36|] do not agree with the more recent 


experiments [40, 41, [18]. Furthermore, the error presented in [36|J appears 
to be too small, especially taking into account that the signal-to-background 
ratio in this experiment is equal to ~ 1/8. It has been mentioned before [77] 
that the half-life value in this work can be ~ 1.5 — 2 times higher because the 
thickness of the dead layer in the Ge(Li) detectors used can be different for 
crystals made from enriched Ge, rather than natural Ge. With no uniformity 
of the external background (and this is the case!), this effect can have an 
appreciable influence on the final result. 

Finally, in calculating the average, only the results of experiments with 
signal-to-b ackg round ratios greater than 1 were used (i.e., the results of Refs. 

ET" 


40 


18]). The weighted average value is: 


Ti /2 = 1.65 


+ 0.14 

- 0.12 


• 10 21 yr. 
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3.3. S2 Se 

There are three independent counting experiments and many geochemical 
measurements (~ 20) for 82 Se. The geochemical data are neither in good 
agreement with each other nor in good agreement with the data from the 
direct measurements. Typically, the accuracy of geochemical measurements 
is at the level of 10% and sometimes even better. Nevertheless, the possibility 
of existing large systematic errors cannot be excluded (see discussion in Ref. 

78]). Thus, to obtain a present half-life value for 82 Se, only the results of 


the direct measurements 
preliminary result of 
result of work 


42 


ll, 42, |43| were used. The result of Ref. 


is the 

hence it has not been used in our analysis. The 
is presented with very asymmetrical errors. To be more 
conservative only the top error in this case is used. As a result, the weighted 
average value is: 

20 , 


T 1/2 = (0.92 ± 0.07) • 10 2U yr. 


3.4. 96 Zr 



a present weighted half-life value for 96 Zr of: 


T 1/2 = (2.3 ± 0.2) • 10 19 yr. 


3.5. 100 Mo 

There are eight positive results from direct experiment^] and one result 
from a geochemical experiment. I do not consider the preliminary result of 
Elliott et al. jloj and instead use their final result j53|, plus I do not use the 
geochemical result (again, see comment in Sec. 3.3). Finally, in calculating 
the average, only the results of experiments with signal-to-background ratios 
greater than 1 were used (i.e., the results of Refs. [sh - firl Q0). In addition, 
I have used the corrected half-life value from Ref. [5l| (see explanation in 
0 ). The following weighted average value for this half-life is then obtained 
as: 

Ti /2 = (7.1 ±0.4) • 10 18 yr. 


1 I do not consider the result of Ref. [80| because of a potentially high background 
contribution that was not excluded in this experiment. 

















3.6. 100 Mo - 100 Ru (T)+; 1130.32 keV) 

The transition to the Of excited state of 100 Ru was detected in seven inde¬ 
pendent experiments. The results are in good agreement, and the weighted 

” 2l| 


average for the half-life using the results from [56], |57|, [591, [60], |20 


is: 


Ty 2 = 6.7 


+ 0.5 

- 0.4 


• 10 20 yr. 


The result from 581 was not used here because we considered the result from 


59] as the final result of the TUNL-1TEP experiment. 


3.1. m Cd 

There are five independent positive results [35], [6j], 63, fhl, [zl that are in 
good agreement with each other when taking into account the correspond¬ 
ing error bars. Again, 1 use here the corrected result for the half-life value 
from Ref. |6j|]. The original half-life value was decreased by ~ 25% 
explanation in 0). The weighted average value is: 


see 


T 1/2 = (2.87 ±0.13) • 10 19 yr. 

3.8. 128 Te and 130 Te 

For a long time, there were only geochemical data for these isotopes. 
Although the half-life ratio for these isotopes has been obtained with good 
accuracy (~ 3%) [66j], the absolute values for Xi/ 2 of each nuclei are different 


from one experiment to the next. One group of authors 65], [81], |82] gives 


T \/2 ~ 0.8 • 10 21 yr for 130 Te and Xi/ 2 ~ 2 • 10 24 yr for 128 Te, whereas another 


group [44], [66| claims Ti/ 2 ~ (2.5 — 2.7) • 10 21 yr and Ti/ 2 ~ 7.7 • 10 24 yr, 
respectively. Furthermore, as a rule, experiments with young samples (~ 
100 million years) give results of the half-life value of 130 Te in the range of 
~ (0.7 — 0.9) • 10 21 yr, while old samples (> 1 billion years) have half-life 


values in the range of ~ (2.5 — 2.7) • 10 21 yr. Recently it was argued that 


short half-lives are more likely to be correct |67l. 1681], Using different young 
mineral results, the half-life values were estimated at (9.0 ± 1.4) • 10 2 ° yr 
and (8.0 ± 1.1) -J.0 20 yr [68f for 130 Te and (2.41 ± 0.39) ■ 10 24 y [67] and 


(2.3 ± 0.3) • 10 24 yr 0 for 128 Te. 


The Erst indication of a positive result for 130 Te in a direct experiment 
More _accurate and reliable value was obtained recently 
The results are in good agreement, and the 


was obtained in 69 


in NEMO-3 experiment 23 


weighted average value for half-life is 


Ty2 = (6.9 ± 1.3) • 10 2 V 
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Now, using very well-known ratio Ti/ 2 ( 130 Te)/ri/ 2 ( 128 Te) = (3.52±0.11)-10 4 
66 ], one can obtain half-life value for 128 Te, 


T y 2 = (2.0 ± 0.3) • 10 M yr. 

I recommend to use these last two results as the best present half-life values 
for 130 Te and 128 Te, respectively. 

3.9. 136 Xe 

The half-life value was recently measured in two independent experiments, 
EXO 28J, 127] and Kamland-Zen 83, [26]. To obtain avera ge value I use most 
precise results from these experiments, obtained in j26], 27] (see Table 1). 
The weighted average value is 

Ti /2 = (2.19 ±0.06) • 10 2 V. 


3.10. 150 Nd 

This half-life value was measured in three independent experiments 


53 


2ft 


71]. The most accurate value was obtained in Ref. [71j]. This value 
is higher than in Ref. [53j and lower than in Ref. [70] (~ 3cr and ~ 2 a 
differences, respectively). Using Eq. (1), and three existing values, one 
obtains Tiy 2 = (8.2 ± 0.5) • 10 18 yr. Taking into account the fact that y 2 > 1 
and S = 1.89 (see Eq. (2)) we then obtain: 

T 1/2 = (8.2 ± 0.9) • 10 18 yr. 

3.11. 150 Nd - 150 Sm (0+; 740.4 keV) 

There are two independent experiments in which 2 v(3f3 decay of 1&0 Nd to 
the 0+ excited state of 150 Sm was observed 72, 24| (the preliminary result 
of Ref. [72] was published in Ref. [84| ). The results are in good agreement. 
The weighted average value is: 


T 1/2 = 1.2 


+ 0.3 

- 0.2 


• 10 20 yr. 


3.12. 238 U 

There is only one positive result but this time from a radiochemical ex¬ 
periment [73|: 

T 1/2 = (2.0 ± 0.6) • 10 21 yr. 
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Table 2: Half-life and nuclear matrix element values for two neutrino double beta decay 
(see Sec. 4). For 130 Ba G 2v value for ECEC transition is taken from [30]. “1 Obtained 
using SSD model. 


Isotope 

7i/ 2 (2^), yr 


1 | 




[G 2v from [29]) 

[G 2v from [32]) 

recommended 





value 

48 Ca 

4.4 + q5 • 10 19 

n n889+ u - uu ' !!4 

U.UOOZ „ 0 0024 

n 0889+0-00^4 
U.UOOZ„q 0024 

0.038 ±0.003 

76 Ge 

1.65^;] 4 • 10 21 

fill OO+0.0043 

U. 11 ZZ„q 0045 

0-1143^:“ 44 

0.113 ±0.006 

82 Se 

(0.92 ±0.07) • 10 2 ° 

o.o826±8:88®; 

0.083ll°:“ 33 

0.083 ±0.004 

96 Zr 

(2.3 ±0.2) • 10 19 

0.0798;°;“ 37 

0.08041S;”! 

0.080 ±0.004 

100 Mo 

(7.1 ±0.4) • 10 18 

0 2065 + ° 0061 
u.zuoo _ 0 0056 

n 18/17+0.0050“) 
U.io4f _ 0 0031 

0.2088^;^ 2 

0.185 ±0.005 

100 Mo- 

6.7loi • 10 20 

0 1 571 +°- 0048 

U.19( 1 —0.0056 

0 161 9 +00050 
u.iury _ 0 0058 


1 °°Ru(0 ] f ) 


n 1 e:iq+0.0047 a ) 

U. 1010 _Q 005 3 


0.151 ±0.005 

116 Cd 

(2.87 ±0.13) • 10 19 

Oil oq+0-0026 

U. 1 Iz-O_o 0024 

0 l04q+°-° 024a) 
u.ruuy _ 0 0023 

O.1139^8:8o25 

0.105 ±0.003 

128 Te 

(2.0 ±0.3) • 10 24 

0.0431±°;“ 37 

0.0483^q;oq34 

0.046 ± 0.006 

130 Te 

(6.9 ± 1.3) • 10 2 ° 

0.0308^;“ 34 

0 0 31 7 + 0 - 0034 

U.UOl 1 -0.0026 

0.031 ±0.004 

136 Xe 

(2.19 ±0.06) • 10 21 

n ni 77 + 0.0003 
U - U1 * '- 0.0002 

0 018c;+ 0 0003 
U.UioO_ Q 0002 

0.0181 ±0.0007 

150 Nd 

(8.2 ±0.9) • 10 18 

0.0579^8:0029 

0 0587 + °' 0034 

U.UOO 1 -0.0030 

0.058 ±0.004 

150 Nd- 

15 O Sm(0+) 

1 . 2±°; 3 • 10 2 ° 

0.0438±°;°° 42 

f) 0450 +a ° 043 
u.uuou-o 0048 

0.044 ±0.005 

238pj 

( 2.0 ± 0 . 6 ) • 10 21 

0 1 8C18+ 0 - 0361 

U.lOOO _ 0 Q227 

0 071 8+ 0 0139 

u.u 1 ro_ 0 o0 88 

n 1 Q+0.09 
u - io -0.07 

i 30 Ba, 

ECEC(2i/) 

~ 10 21 

~ 0.26 [30] 


~ 0.26 


3.13. 130 Ba (ECEC) 

For 130 Ba positive results were obtained in geochemical measurements 
only. In geochemical experiments it is not possible to recognize the different 
modes. But I believe that exactly ECEC(2z/) process was detected because 


other modes are strongly suppressed (see, for example, estimations in 11 


85 


74 , in which 


86 ]). First positive result for 130 Ba was mentioned in Ref. 
experimental data from Ref. [87] were analyzed. In this paper positive result 
was obtained for one sample of barite (T 1/2 = 2 . 11 q; 8 - 10 21 Y r )> but for second 
sample only limit was established (T 1/2 > 4 • 10 21 yr). Then more accurate 


half-life values, (2.2 ± 0.5) • 10 21 yr [75] and (0.60 ± 0.11) • 10 21 yr [25 


were 
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obtained. However, the results are in strong disagreement. One can not use 
usual average procedure in this case. One just can conclude that half-life 
of 130 Ba is ~ 10 21 yr. To obtain more precise and correct half-life value for 
130 Ba new measurements are needed. 


4. NME values for two neutrino double beta decay 


A summary of the half-life values are presented in Table 2 (2-nd column). 
From the measured half-life one can extract the experimental nuclear matrix 


element using the relation 29 


Tyl = G 2u 


9 A 


m P c 


M- 


2u) 


(3) 


where Ti/ 2 is the half-life value in [yr], G 2v is the phase space factor in 
[yr -1 ], gA is the dimensionless axial vector coupling constant and (;m e c 2 * ■ M 2u ) 
is the dimensionless nuclear matrix element. It is necessary to take into 
account that there are various indications that in nuclear medium the matrix 
elements of the axial-vector operator are reduced in comparison with their 
free nucleon values. This quenching is often described as a reduction of 
the coupling constant gA from its free nucleon value of q& = 1.2701 76 
to the value of gA ~ 0.35 — 1.0 (see discussions in [i, 0, 0, 0]). So, follow 


the Ref. |29| it is better to have a deal with so-called ’’effective” NME, 
,2 - ' I And this value has been calculated for all 


I M 2u S I = 9A' I ( m e c ' Z ■ M 2 u 
mentioned above isotopes. 

The results of these calculations are presented in Table 2 (3-d and 4- 


29 


th columns). To do the calculations I used the G 2u values from Ref. 

31, 32 H, respectively (see Table 3). For 130 Ba G 2v value for ECEC 


and 

transition was taken from |30j]. These recent calculations pretend to be most 


d, |3ll, 1320) 


reliable and correct by this moment (see discussions in [29], 

Results of these calculations are in quite good agreement (~ 1-7%) with two 
exceptions, for 128 Te (~ 20%) and 238 U. For 238 U two absolutely different 
values 14.57 • 10 -21 yr -1 (29j and 98.51 • 10 -21 yr -1 [32]) were obtained. It 
is clear that calculations for 238 U have to be checked. For 100 Mo, 100 Mo- 
100 R/u(Oj~) and 116 Cd I used G 2u calculated in Ref. |29| for SSD mechanism, 


in addition. Corresponding values for 




are presented in Table 2 


2 Ref. [32] was published as up-date of Ref. l31]. And, finally, I used in this work results 

of calculations from Ref. 1321. 
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and these are most correct values for these isotopes. So-called recommended 
values for | M^l | are presented in Table 2 (5-th column) too. These values 
were obtained as an average of two values, specified in columns 3 and 4. The 
error of recommended values is chosen so that to cover all range of values 
from columns 3 and 4 (taking into account corresponding errors). For 100 Mo, 
100 Mo- 100 Ru(0 1 1 ') and 116 Cd I recommend to use values obtained using Giv 
for SSD mechanism. 

For the majority of isotopes we now have | Mjy | with an accuracy of 
~ 3 — 8%. For 128 Te and 130 Te it is ~ 13% and for 150 Nd- 150 Sm(Oj h ) it is 
~ 16%. The most unsatisfactory situation is for 238 U (~ 70%) and 130 Ba 
(~ 50%). For 238 U main uncertainty is connected with accuracy of G 2i , and 
for 130 Ba with accuracy of experimental data for the half-life. 

In a few recent publications 0,1,0 attempts to reproduce NME values 
for two-neutrino double beta decay within various models were realized. The 
conclusion was that renormalization (quenching) of qa is needed to reproduce 
the experimental data. So within Interacting Shell Model (ISM) approach 0] 
NMEs for 48 Ca, 76 Ge, 82 Se, 128 Te, 130 Te and 136 Xe were calculated, using the 
gA ~ 0.57 — 0.94 (these values were obtained from data for a single beta decay 
or charge exchange reactions). As a result it was succeeded to obtain rather 
good agreement between calculations and experimental data (nevertheless, 
NME calculated values for 82 Se, 128 Te, 130 Te and 136 Xe exceed experimental 
data on r\j 20 — 30%). In Ref. 0] within QRPA model calculated values of 
NMEs for 100 Mo, 116 Cd and 128 Te were adjusted to experimental values at 
the expense of a choice of the corresponding gA values (~ 0.4 — 0.75). The 
same procedure was executed within IBM-2 model for many nuclei [hj. It 
was shown that for exact reproduction of experimental data the gA has to 
be ~ 0.35 — 0.71. The question of whether or not the quenching of gA is the 
same in 2 u/3/3 as in 0 u/3/3 decay is the subject of debate, but it is clear that 
this question has to be carefully investigated because changes in gA leads to 
changes in sensitivity to effective Majorana neutrino mass (m u ) in double 
beta decay experiments. 

I would like to note that in all these cases 0,0,0 when comparing with 
experimental data the recommended 7\/ 2 ( 2u) values from our previous work 
|l7j | were used. In the present work more precise experimental values for 
Ti/2(2v) and NME for many nuclei are obtained and, I believe, that will help 
with a solution of the gA problem in the future. 
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Table 3: Phase-space factors from Ref.[29|, @ and do|. “) Obtained using SSD model. 


Isotope 

GW(io- 2 V 1 ) [22] 

G 2u (10- 21 yr~ 1 ) [32] 

48 Ca 

15550 

15536 

76 Ge 

48.17 

46.47 

82 Se 

1596 

1573 

96 Zr 

6816 

6744 

100 Mo 

3308 

4134“^ 

3231 

io °Mo- 100 Ru(0+) 

60.55 

65.18“) 

57.08 

116 Cd 

2764 

3176“) 

2688 

128 Te 

0.2688 

0.2149 

130 Te 

1529 

1442 

136 Xe 

1433 

1332 

150 Nd 

36430 

35397 

150 Nd -150 Sm (Q+) 

4329 

4116 

238pj 

14.57 

98.51 

130 Ba, ECEC(2z/) 

15000 [30] 



5. Conclusion 


In summary, all positive 2z/ / 9/3-decay results were analyzed, and average 
values for half-lives were calculated. For the cases of 12S Te and 130 Ba, the so- 
called recommended values have been proposed. Using these half-life values, 
I | for two neutrino double beta decay were obtained. Finally, previous 

results from Ref. 17] were successfully up-dated and new results for 136 Xe 


were added. A summary is collected in Table 2. I strongly recommend the 
use of these values as the most reliable presently. 

Notice that the accurate half-life values for 2 u(3j3 decay could be used to 
adjust the most relevant parameter of the quasiparticle random-phase ap¬ 
proximation (QRPA) model, the strength of the particle-particle interaction 
g pp . In addition effective qa value could be established for 2/3 decay. It will 
give the possibility to improve the quality of NME calculations for neutrino¬ 
less double beta decay and, finally, to improve the quality of neutrino mass 
(m„) estimations. 
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